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Quantum Monte Carlo methods are sophisticated numerical techniques for simulating interacting 
quantum systems. In some cases, however, they suffer from the notorious "sign problem" and 
become too inefficient to be useful. A recent publication [J. Wojtkiewicz, Phys. Rev. B 75, 174421 
(2007)] claims to have solved the sign problem for a certain class of frustrated quantum spin systems 
through the use of a bipartite valence bond basis. We show in this Comment that the apparent 
positivity of the path integral is due to a misconception about the resolution of the identity operator 
in this basis, and that consequently the sign problem remains a severe obstacle for the simulation 
of frustrated quantum magnets. 



Frustrated quantum antiferromagnets are of great in- 
terest, from both experimental and theoretical points of 
view. A theoretical analytical framework for understand- 
ing frustrated quantum magnets is still lacking, and nu- 
merical simulations also turn out to be very difficult. In 
particular. Quantum Monte Carlo (QMC) methods suf- 
fer from the "sign problem," rendering QMC simulations 
of frustrated magnets completely inefficient. 

In a recent publication, ^ J. Wojtkiewicz claims to have 
solved the sign problem for a certain class of frustrated 
antiferromagnets, including the important cases of the 
J1-J2 model and the pyrochlore lattice. He proposes us- 
ing a valence bond (VB) basis, rather than the conven- 
tional Sz basis, to formulate the path integral used in the 
QMC framework. In this basis, the matrix elements in- 
volved are strictly positive, and the sign problem seems 
to disappear. In this Comment, we point out that the 
conventional route to constructing the path integral (via 
the Trotter- Suzuki formula) is not suited to the bipartite 
VB basis used by Wojtkiewicz, and that consequently the 
method proposed in Ref. Jj cannot work. 

Let us first summarize the approach taken in Ref. [H. 
The starting point is to group terms in the Hamiltonian 
H = J2p Hp according to the plaquette P in which they 
act. The Trotter-Suzuki formula Isee Eq. (1) in Ref. Ill, 
is used to decouple the potentially noncommuting terms 
appearing in the partition function: 
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Here, t — (3/n is an infinitesimal imaginary-time step. 
In the usual way, a representation of unity is introduced 
between each exponential factor, leading to 
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where the notation implies \ap^n+i) = \oip.i)- This 
amounts to introducing a new set of states at each pla- 
quette and at each time slice. 



Implicit in Eq. ^ is the existence of a simple resolution 
of the identity operator. 
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in the basis {a} used to perform the simulations. Taking 
this for granted, Ref. Hi calculates all possible matrix ele- 
ments {a\e~'^^^\a') (with la) and \a') in the bipartite VB 
basis) and shows that they are all positive. Wojtkiewicz 
then concludes that a QMC process without sign problem 
can be devised with this decomposition. 

The essence of this Comment is simple: we point out 
that Eq. ^ does not hold in the bipartite VB basis con- 
sidered for the calculations. The correct resolution of 
identity reintroduces the sign problem in the Trotter- 
Suzuki decomposition of the partition function and ren- 
ders the scheme proposed in Ref. 1 useless. 

In the bipartite valence bond (BVB) basis, the lattice 
is decomposed into two sets (A and B) of sites, and the 
two spins forming a singlet necessarily belong to differ- 
ent sets: [a, 6] = (ITaib) - |iaTfc))/\/2- The bipartite 
VB basis is overcomplete, albeit smaller than the full va- 
lence bond (FVB) basis, where there are no restrictions 
on the constituent spins of the singlet. Since VB bases 
are nonorthogonal, the most general form of the identity 
operator is 
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This means that one should also consider the matrix el- 
ements of la,a' at the QMC sampling level, and in par- 
ticular ensure that they are all positive in order to avoid 
the sign problem. 

The massive overcompleteness of the basis offers mul- 
tiple different solutions for la.a'- One of the possible 
choices in the FVB basis is particularly simple. Indeed, 
for a system of 2N spins, la.a' can take the diagonal form 
[^^^ ^'l- (27) in Ref. 2]. But, such a simple 
resolution is not possible if the basis consists of AB bonds 
only3 (BVB basis). 
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To see this explicitly, consider the four-site (A'' = 2) 
system, which admits three VB states — two bipartite, 
[1,2][3,4] and [1, 4] [3, 2], and one nonbipartite, [1,3][2,4]. 
The representation of the identity operator has only pos- 
itive, diagonal elements in the full VB basis but negative, 
offdiagonal ones in the bipartite basis: 

/I 0\ 

(1.,o')fvb = 3 ^ (5) 

Let us comment on the general N case: the matrix 
1q,c(' is the pseudo-inverse of the positive semidefinite 
overlap matrix Oa,a' = {o:\a') (when only AB bonds are 
allowed). Since Oa.a' has only positive elements, the 
positivity of all off-diagonal elements of la, a' would re- 
quire that all diagonal elements be strictly negative. This 
is in condradiction with the positivity of the traces of 
both Oct, a' and la,a'- Therefore, in the basis of bipartite 
bonds, la^a' must contain at least one strictly negative, 
off-diagonal element. Indeed, one can establish that la. a' 
is dense in such elements. 

Alternatively, one could imagine formulating the QMC 
using the full VB basis, which has a simpler resolution of 
identity. Unfortunately, the sign problem then appears 
in the Hamiltonian matrix elements. In short, we are left 
with two ways to organize the path integral, one involving 
a sum over a VB basis restricted to AB bonds only, 
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and one involving an unrestricted sum, 

FVB 
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In neither case is the sign eliminated. In the bipar- 
tite VB basis, the matrix elements {a\e~'^^^\a') are 
strictly positive, but the identity matrix elements la. a' 
are not. In the full VB basis, the situation is reversed: 
lct,Q' ~ Sa,a' > 0, but {a\e~'^^^\a') is no longer strictly 
positive. Thus, the choice of basis amounts to nothing 
more than shuffling the sign problem from one location 
to another. With these considerations in mind, it is clear 
that the Trotter-Suzuki decomposition cannot be per- 
formed as envisioned in Ref. [l|, and the QMC scheme 
proposed therein breaks down. 

In conclusion, the VB basis can be very useful in QMC 
simulations, but a workable algorithm cannot be orga- 
nized around inserting an overcomplete set of states be- 
tween time slices in a Trottcr-Suzuki scheme. A better 
approach is to project down from a trial state using the 
power method, relying on the property that Heisenberg 
spin interactions induce simple rearrangements of valence 
bond configurations.^^ In this way, spin models with non- 
frustrating interactions of arbitrary order can be treated. 
Grouping of terms (along the lines of the plaquette de- 
composition of Ref. 1) combined with a fixed- node ap- 
proximation^ might partially alleviate the sign problem 
in models with frustrating interactions and lead to more 
efficient QMC sampling,^ but we do not expect a simple 
full solution to the sign problem (such as the one claimed 
in Ref. 11) to emerge from this approach. 

Finally, it is worth noting that the overlap calculations 
reported in Ref. 1, are certainly correct and may yet prove 
useful in the evaluation of variational valence bond wave- 
functions, following the method of Liang, Dougot, and 
Anderson^ and its refinementsi^*^ 
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